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Heavy quark bound states are used as significative probes of the QCD vacuum and 
the mechanism of confinement. 



1 Strong interaction and topological nontrivial configurations 

The focus of this school is on the "Topology of strongly correlated systems" 
in different areas ranging from QCD to condensed matter systems. In all 
these fields, investigations were presented about nontrivial topological configu- 
rations. In QCD such configurations are believed to be related to thej^nixivial 
structure of the QCD vacuum and to the confinementa^echanism EJIj'qou on 
one hand and to the breaking of the chiral symmetry uuu on the other hand. 
However, the QCD dynamics is extremely complicate and in order to test our 
understanding of the topological configurations, we need a systematic and un- 
der control way of parameterize the low energy physics in a frame that can still 
be simply related to the real experiments and/or to the lattice experimental. 
It is the aim of my talk to show how this is possible and in which physical 
situations. 

2 Heavy Quark Bound States 

Heavy-quark bound systems are an ideal playground for theoretical ideas about 
strong QCD. In fact, while light quarks are highly "nonperturbative objects", 
strongly interacting with the vacuum and with a mass almost entirely domi- 
nated by chiral symmetry breaking (topological nontrivial) effects, heavy quarks 
behave as "external sources" and thus as natural probes of the QCD vacuum. 
Being heavy at least the quark mass scale m can be treated perturbatively: 
Tfi 3> Aqcd- Aqcd has here always to be understood as the scale at which 
nonperturbative effects become dominant. Furthermore, since the characteris- 
tic difference AE between the energy levels is AE <C m inside such systems, 
they are nonrelativistic and can be described at leading order by the appropri- 
ate Schrodinger equation. This gives origin, as usual in nonrelativistic bound 
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Figure 1: Tqq vs Cornell potential (from Ref. s ). 

state, to several dynamical scales below m: the relative momentum p ~ mv 
(the inverse of which gives the characteristic spatial size r of the system), the 
energy E ~ mv 2 (the inverse of which gives the characteristic time T of the 
system). Here, v is the velocity of the heavy quark in the bound state. In 
QCD, v is in general a function of the whole dynamics, perturbative and non- 
perturbative, but still remains a small parameter, v -C 1, and thus the scales 
mentioned above turn out to be well separated and hierarchically ordered. It 
is well known that a technical (and quite hard) problem in bound state calcu- 
lation is originated by the fact that such scales get entangled EJ. Thus, even 
in QED it turns out to be very useful to disentangle such scales. In QCD 
this becomes a deeper and more conceptual issue, since we need to separate 
as much as possible the perturbative contributions ('hard' physics at the high 
scale) from the nonperturbative effects ('soft or ultrasoft' physics at the low 
scale). It is on these last contributions that we will eventually be in the position 
to use our knowledge about topologically nontrivial configurations. 

I will show here that in quarkonium physics the existence of different dy- 
namical scales besides Aqcd 

• makes the calculations more complicate from a technical point of view; 

• however, it makes also possible to test the non-perturbative nature of the 
QCD vacuum at different levels of deepness. 

In fact, heavy quark bound states provide a full set of probes that may ex- 
plore different characteristic distances, from the small to the large distance 
regime, and different dynamical situations, thus providing insights into pertur- 
bative/non-perturbative factorization mechanisms. In Fig. 1 we show a phe- 
nomenological determination of the characteristic radius r of heavy quarkonia 
against the phenomenological QQ potential. It is apparent that heavy quarko- 
nia display a pattern of characteristic radii that extend from the perturbative 
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Figure 2: The more extended the physical object, the more sensitive to non-perturbative 
physics. 

region r <C Aqcd to the transition and the confinement region r ^> Aqcd- 
For ground state quarkonia (e.g. T(15)) r <C 1/Aqcd: the scale 1/r is per- 
turbative (this also implies that the binding potential is perturbative). For 
most quarkonia r ~ 1/Aqcd: the scale 1/r is nonperturbative. In order to 
be able to perform a systematic study, we need a clean and under control 
approach that allows us to factorize and resum the perturbative physics at 
the appropriate scale, leaving us with the appropriate parameterization of the 
nonperturbative physics at the low energy scale. This goal is achieved by con- 
structing an effective field theory (EFT) . The existence of a hierarchy of energy 
scales in quarkonium systems allows the construction of EFT with less and less 
dynamical degrees of freedom but completely equivalent to QCD. This leads 
ultimately to a field theory derived quantum mechanical description of these 
systems. We call potential nonrelativistic QCD (pNRQCD) the corresponding 
EFT. pNRQCD provides an unambiguous power counting in v which deter- 
mines w^chxiperators are relevant at a given order of the expansion in v. We 
refer toB-BEl for details about pNRQCD. Here, we concentrate on the non- 
perturbative contributions: depending on the extension of the 'probe' systems 
such contributions enter parameterized in a completely different way. This 
is precisely the interesting point that may allow us to understand something 
about the QCD vacuum structure. In QCD we can never get completely rid of 
the nonperturbative effects. However, I will show that a pattern like the one 
represented in Fig 2 is realized in the heavy quark bound state nonperturba- 
tive dynamics. If the spatial and temporal characteristic scales of the bound 
system are perturbative (i.e. if mv 3> mv 2 3> Aqcd)> then nonperturbative 
contributions arise only in the form of local gluon condensates which contribute 
to the energy levels. When the characteristic temporal scale starts to be dom- 
inated by Aqcd {mv 2 ~ Aqcd), then the nonperturbative contributions are 
encoded in nonlocal (time-dependent) gluon condensates. As soon as 1/Aqcd 
falls between the spatial and the temporal scale (mv > Aqcd > mv 2 ), the 
nonperturbative nonlocal condensates enter the definition of the potential and 
give origin to short distance (r < 1/Aqcd) nonoperturbative corrections. At 
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the end, when Aqcd is comparable to the spatial size (rnv ~ Aqcd), the non- 
perturbative information is carried by bidimensional extended objects called 
Wilson loops. The potential is given in terms of them. Quite interestingly, the 
Wilson loop turns out to be the order parameter of confinement (in gluody- 
namics) and allows us to put in direct connection topological configurations 
and the heavy quark phenomenology. 

As we see from Fig. 2, there are sufficient physical systems in nature to 
test all these dynamical situations. 

3 Nonrelativistic quark bound states with a small characteristic 
radius: Condensates 

For heavy quark bound systems for which the following hierarchy m 3> mv 3> 
I7iv > Aqcd holds, we can sequentially integrate out in perturbation theory 
both the hard scale m and the soft scale mv. Then v ~ a s . We denote by 
R = (xi + xa)/2 the center of mass of the QQ system and by r = xi — X2 
the relative distance. Then at the scale of the 'matching' // (mv 3> fi' 3> 
tow 2 , Aqcd) we have as still dynamical degrees of freedom: QQ states with 
energy of order of the next relevant scale, 0(Kqcd, mv 2 ), momentum of order 
0(mv), plus ultrasoft gluons A M (R, t) with energy and momentum of order 
0(Aqcd> mv ). Therefore only very low energy gluons remain in the EFT as 
dynamical degrees of freedom: they are multipole expanded. We find convenient 
to decompose the QQ state into a singlet 5(R, r, t) and an octet 0(R, r, t) with 
respect to color transformations. The pNRQCD Lagxangian is given at next 
to leading order (NLO) in the multipole expansion bji3: 

£ = Tr{ S t (^ _E!_T4(r)-£^vWot ( iDo - ^ - V (r) - £ ^) o) 

+gV A (r)Ti {O f r • ES + S f r • EO} + g^p-Tr {0*r • EO + O^Or ■ E} - ^F^F^ ua . (1) 

All the gauge fields in Eq. (|l|) are evaluated in R and t. In particular 
E = E(R, t) and iD O = idoO — g[Aa(R,, t), O]. The Vj are potential matching 
coefficients: they are in the EFT the correspondent to the usual notion of the 
potential. We call V s and V a the singlet and octet static matching potentials 
respectively. The singlet sector of the Lagrangian gives rise to equations of mo- 
tion of the Schrodinger type, while the last line of (|l|) contains retardation (or 
nonpotential) effects that start at the NLO in the multipole expansion. At this 
order the nonpotential effects come from the singlet-octet and octet-octet inter- 
action mediated by an ultrasoft chromoelectric field. Recalling that r ~ 1/ mv 
and that the operators count like the next relevant scale, 0(mv 2 , Aqcd), to the 
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power of the dimension, it follows that each term in the pNRQCD Lagrangian 
has a definite power counting. In the EFT language the potential is defined 
upon the integration of all the scales up to the ultrasoft scale mv 2 and thus if 
Aqcd <- fnv 2 the potential is given by a pure perturbative expansion in a s at 
all orders. The EFT gives us with (jl]) a precise prescription to calculate both 
potential and nonpotential contributions to the energy levels and to count their 
relevance in power of a s . Nonpotential effects start only at order ot 5 s ln/i'. At 
this order (three loops) the static singlet potential starts_Lo become sensitive 
to the infrared physics via a logarithmic dependence on /iEU. Such dependence 
is re-absorbed in a physical observable by non-potential contributions^. Non- 
perturbative effects are purely of nonpotential nature and the leading ones are 
carried by the chromoelectric field which mediates the singlet-octet interaction. 
From (|j) thfLfluarkonium energy levels have been calculated up to a|? In // and 
are given b]£L3 



to to 

2 roo 

"W^/o d *( n ^i rei4(E " _Ho)r l n ^)< EtI (*)^(*' )«6 jE6 (0)>(-"') (2) 

C 2 a 2 2 (0) 

where E n — —to £ £ , H a = E- + Vo and \n, I) are the Coulomb wave func- 
tions. V s is taken at three loop leading log and the other singlet potentials in 
the 1 /to expansion at the appropriate order in a s for the requested accuracy. If 
Aqcd <C mv 2 , also the scale mv 2 can be integrated out perturbatively from the 

above formula and, noticing that: m j^c7oTY = 3N~Tf ( n i I r e -h 
the non-perturbatijze contributions reduce simply to the well-known Voloshin 



Leutwyler formula 



. j e n n^ 2 (F 2 (0)) 



(e n is a (known) number of order 1) and thus the only nonperturbative con- 
tribution is given in this case at leading order by the local gluon conden- 
sate (F 2 (0)) = {{a s /ir)F^ v (0)F a ^(0)). This is a well known nonperturbative 
parameter that can be calculated also inside QCD vacuum models Fur- 
ther nonperturbative corrections are still of nonpotential type and are sup- 
pressed. For quarkonium of a typical size smaller than 1/Aqcd, the most 
relevant operator for the nonperturbative dynamics is the bilocal gluon con- 
densate (E a (t)4>(t, 0)^E fc (0)), which belongs to the class of th non-local non- 
perturbative gluon condensate (F^ l/ (x)<p(x, 0)* b *F* „(G)}- Let me open a small 
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parenthesis by summarizing our knowledge of such operator and what model 
independent information can be obtained by the EFT 

Different parametrizations have been proposedllsESEJ'Ea for the nonlocal 
gluon condensate. Because of its Lorentz structure, the correlator is in general 
described by two form factors. A convenient choice of these consists in the 
chromoelectric and chromomagnetic correlators: 

The strength of the correlators is of the order of the gluon condensate. In the 
long range (x 2 — > oo) they fall off exponentially (in the Euclidean spacel_with 

0me typical correlation lengths that can be measured on the lattice EJjO. In 
two different lengths have been obtained for the chromoelectric (T^) and 
the chromomagnetic correlators [Tg 3 ): 

Tg + Tg 3 ~ 0.1-0.2 fm (quenched). (4) 

A recent sum-rule calculation @ gives < Tg . The sum rule turns out not 
to be stable for the chromoelectric correlator, while for the chromomagnetic 
correlation length it gives 

T g B = 0.1ll8;8| fm (quenched). (5) 

The interesting thing is that these and , which control the form of the 
nonperturbative correction entering the energy levels, have a precise physical 
interpretation. Their inverses correspond to the masses of the lowest-lying 
vector and pseudovector gluelumps^, respectively. This can be explicitly seen 
in the short-range liraitp-r = |x| — > 0, where the hybrids operators can be 
explicitly constructed The suitable effective field theory is pNRQCD in 
the static limit. Gluelump operators are of the type Tr{OH}, where O = O a T a 
corresponds to a quark-antiquark state in the adjoint representation (the octet) 
and H = H a T a is a gluonic operator. By matching the QCD static hybrid 
operators into pNRQCD, we get the static energies (also called potentials) of 
the hybrids, which reduce to the gluelumps in the limit r — > 0. At leading 
order in the multipole expansion, they read 



V H (r) = V {r) + - w + 0(r 2 ); (H a (t)<P(t, )^ H b (0)}™-^ ~ k" 1 ^ 

9 



(G) 



b The gluelump is a state formed by a QQ pair in the adjoint representation plus a gluon, all 
in the same spatial point. They are measured on the lattice and their masses are, roughly 
speaking, the limit of the hybrids static energies when the typical QQ distance r goes to 
zero. 
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Since hybrids are classified in QCD according to the representations of -Doo.h, 
while in pNRQCD, where we have integrated out the length r, their classifi- 
cation is done according to the representations of 0(3) x C, the static hybrid 
shontpiige spectrum is expected to be more degenerate thanJdie long-range 
oneEj'El. Tliii lattice measure of the hybrid potentials done in til confirms this 
feature. In EU it has been shown that the quantum numbers attribution of 
pNRQCD to the short-range operators, and the expected 0(3) x C symmetry 
of the effective field theory match the lattice measurements. By using only E 
and B fields and keeping only the lowest-dimensional representation we may 
identify the operator H for the short-range hybrids called S+ (and IT g ) with 
r • E (and r x E) and the operator H for the short-range hybrids called E~ 
(and n u ) with r • B (and r x B). Hence, the corresponding static energies for 
small r are 

9 a 

The lattice measure of0 shows that, in the short range, VL+ < „ (r) > V s - n (r) . 

r~l .g • y " ' u 

This supports the sum-rule predictionllj that the pseudovector hybrid lies lower 
than the vector one, i.e. Tg < Tj? . 

Summing up, the message is that, within the EFT, we can obtain quite 
a number of model-independent information on the nonperturbative contri- 
butions to the energy levels which are carried by these nonlocal condensates 
and we can establish interesting relations between these objects and physical 
'nonperturbative' entities like the hybrids. 

Let us now briefly consider the situation when mv > Aqcd > rnv 2 . The 
infrared sensitivity of the quark-antiquark static potential at three loops liil 
signals that it may become sensitive to non-perturbative effects if the next 
relevant scale after mv is Aqcd- Indeed, in the situation mv ^> Aqcd 3> 
tow 2 , the leading non-perturbative contribution (in a s and in the multipole 
expansion) to the static potential reads EJ 

2 2 poo 

T/ (r) nonpert = -i^T/- / dte- itC * a °K 2r HE a (t)<l>(t,0)^E b (0))( t i'). (7) 

J'c o Jo 

This term has to be summed to V Q pcTt and it explicitly cancels the dependence 
of the perturbative static potential on the infrared scale \i! . It is interesting to 
note that the leading contribution in the Aqcd / 'mv 2 expansion of y non -p° rt 
(obtained by putting the exponential equal to 1) cancels the order AQ CD r 2 
renormalon that affects the static potential (the leading-order renormalon, of 
order Aqcd, cancels against the pole mass). Therefore, also in the renormalon 
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language, the above operator is the relevant non-perturbative contribution to 
the static potential in the considered kinematic situation. 

The nonperturbative terms in the potential are in this case organized in 
powers of a s /(rAQCD), starting with a quadratic term, r 2 . The actual form 
of these short range nonperturbative corrections to the potential depends on 
the actual form of the nonlocal correlato£il. It is therefore very interesting to 
calculate this object inside models of the QCD vacuum, or with topological 
configurations on the lattice, also in relation to the recent claim of violation of 
the OPE expansion 'carried' by a nonperturbative short distance 'string' term 
in the static singlet potential. 

4 Nonrelativistic quark bound states with a large characteristic 
radius: Wilson loops 

Fig. 1 shows that most of the quarkonia states lie in a region where the inverse 
of the size of the system is close to the scale Aqcd ■ In this situation the poten- 
tial can no longer be expressed as an expansion in a s . The non-perturbative 
dynamics is already switched on at the binding (potential) scale and it is con- 
tained in more extended objects than (local or non-local) gluon condensates: 
Wilson loops and chromoelectric and chromomagnetic field insertions on_the 
Wilson loops. Such operators can be eventually calculated on the lattices or 
in QCD vacuum modelsE5EJ. InE3 the matching of QCD to pNRQCD has been 
performed at order 1/m 2 in the general situation Aqcd < fnv. This has been 
proved to be equivalent to compute the heavy quarkonium potential at order 
1/m 2 . More precisely, a pure potential picture emerges at the leading order 
in the ultrasoft expansion under the condition that all the gluonic excitations 
(hybrids) have a gap of 0(AQCD)|-2fith respect to the singlet. Such situation 
is confirmed by lattice simulations^! Then if we switch off the light fermions 
(pure gluodynamics), only the singlet survives and pNRQCD reduces to a 
pure two-particle NR quantum-mechanical system. Therefore, the situation 
assumed by all potential models, may be nowjgorously derived under a spe- 
cific set of circumstances (for more details sec 2 ]). In this situation, pNRQCD 
involves only a bilinear colour singlet field, iS(xi,X2,t), 



where h s is the Hamiltonian of the singlet, pi = — ?'V Xl and P2 = — «V X2 . It 
has the following expansion up to order 1/m 2 




(8) 




+ V® + 



y(l,0) y(0,l) y(2,0) 00,2) 
+ + — + — 

mi mi m\ m% 



+ 



m\mi 



(9) 
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Higher pccLec effects in the 1/m expansion as well as extra ultrasoft degrees of 
frcedonOH'tJ, such as hybrids and pions can be systematically included and 
may eventually affect the leading potential picture (like in the perturbative 
regime ultrasoft gluonsEJtS). We shall use the following notations: (. . .) will 
stand for the average over the Yang-Mills action, Wn for the rectangular static 
Wilson loop of dimensions r xT and ((...)) = {■ ■ ■ W a ) / (W n ) ■ We define 
((Oi(U r ]...O n (t n ))) c as the connected Wilson loop. We get in terms of Wilson 
loops t3 

V<-°\r)= lim Un(Wn), (10) 

1 — >OC 1 

V^°\r) = ~ lim / T dti(( fl E 1 (t,x 1 )- ff E 1 (0,x 1 )» c = V(°. 1 )(r). (11) 

All the obtained expression for the potentials are manifestly gauge invariant 
and are also correct at any power in a s in the perturbative regime. For the 
terms of order 1 /m 2 we have the general decomposition 

V {2fi) = \ {P?, ^(r)} + ^J^L 2 + V^\r) + V^\r)U ■ S 1; (12) 

where Lj = r x p 3 . Similar definitions hold for V^- ' 2 ^ and V^ 1,1 ). All the poten- 
tials (i.e. the functions of r) in ^(2,0)^(0,2)^(1,1) 

are obtained factorized in 

two contributions. The first one contains the hard physics at the scalc-jp-.which 
is calculated as a series expansion in a s at the appropriate hard scalc3'tPJ. The 
second contribution at the scale 1/r , enjoys a closed expression in terms of 
average values of two or more chromoelectric and chromomagnetic insertions 
in the presence of the Wilson loopEi, similarly to (pX|) , (p"T|) . The explicit ex- 
pressions can be found inl23. Therefore, the calculation of the dynamics, i.e. 
the calculation of the potential, is in this way reduced to a calculation in pure 
gluodynamics at the lowest scale. On one hand, such objects are very conve- 
niently evaluated on the lattice since they are pure glue objects and contain 
only one scale. On the other hand, such objects can be very well evaluated 
iisipg nontrivial topological configurations inside a model of the QCD vacuum 
□3. We emphasize that, since the potential we get here is a well defined 
quantity, derived from QCD via a systematic and unambiguous procedure, 
and complete up to order 1/m 2 , it is not affected by the usual ambiguities 
(ordering, retardation corrections, etc.), which affect all potential models and 
all phenomenological reductions of Bethe-Salpeter kernels For the same 
reason the above results are the relevant ones for the study of the properties 
of the QCD vacuum in the presence of heavy sources. The average value of 
the field insertions in the presence of the Wilson loop contain all the relevant 
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information about the heavy quark dynamics. This is much more detailed and 
quantitative with respect to the generic request that topological configuration 
should originate an area law behaviour into the Wilson loop. If we calculate 
these field insertions on the Wilson loop using some dominant topological con- 
figurations we have then a simple procedure to relate the result, on one hand 
to standard lattice simulations of the same objects and, on the other hand, 
directly to the spectrum. This last feature descends from the fact that within 
the EFT we are able to fully take into account the perturbative contributions 
and we have a clear power counting that selects the terms contributing at a 
given order in the v expansion. 

5 Topology and heavy quark bound states 

The outlined study is rigorous and allows a systematic disentanglement of 
the high from the low energy scales of the heavy quarkonium system under 
study. In this specific sense also perturbative and nonperturbative effects turn 
out to be disentangled. Thus the nonrelativistic dynamics of heavy quark 
bound states is surely an adequate benchmark to test ideas and models of 
QCD topology and QCD vacuum structure. 

As I have shown, when mv ~ Aqcd, then the whole nonperturbative QQ 
interaction at order 1 /m 2 turns out to be expressed in terms of average values 
of Wilson loops and field insertions into the Wilson loop. We-know that in the 
confined phase the Wilson loop is dominated by an area lawJ which is related 
to the formation of an interquark confining flux tubfll. Such feature has been 
obtained using different topological configurations: gluon Abelian projected 
fieldfl, monopolefH, center vorticell. In all these cases an approximated area 
law is obtained and it is debated what are the relevant topological configura- 
tions that really drive the confinement dynamics. Many lattice investigations 
have been performed and have ended up with the notion that in a way all 
these topological confimpations are relatedES. Several QCD vacuum models 
exist in the literaturaJj'ErEjt] and all of them obtain an area law for the Wilson 
loop. However, there are quite a number of differences between different vac- 
uum models as well as between different topological configurations, but such 
differences are not apparent at the level of the static Wilson loop. In this 
case, typically we need only one parameter, the string tension a, to encode the 
information of a constant energy density in the confining flux tube. However, 
there are many ways in which this can be realized in relation to the actual 
profile of the flux tube. It is well known for example that while the QQ flux 
tube measured on the lattice in QCD is quite fat, the one measured on the 
lattice in the Abelian projection is very thirEI. 
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On the other hand in order to describe the nonperturbative heavy quark 
dynamics at order 1 /m 2 we need not only the Wilson loop area law behaviour 
but also the behaviour of the chromoelectric and chromomagnetic field inser- 
tion into the Wilson loop. These last contain information on the profile of the 
flux tube and turn out to be much_more sensitive to the QCD vacuum models 
or the topological configurationtra ^ . Thus, it is much more interesting to 
study the full non relativistic dynamics of heavy quarks in place of the much 
more qualitative statement about finding out confining configurations into the 
Wilson loop. In particular, it would be interesting to have a lattice_calcula- 
tions of (|l0|), (|i~L|), ( p^ ) and the other 1/m 2 poteptials contained incd, using 
monopoles configurations versus center vorticejaO. In conclusion, I would like 
to invite the people working in the field to calculate local gluon condensates, 
nonlocal gluon condensates, Wilson loops and field insertions into the Wilson 
loop in their topological models and to relate them to the phenomenological 
data using the unified and systematic frame of the EFT that I presented here. 
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